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Abstract.  Behaviour  of  equatorial  spread  F  (ESF)  fossil 
plumes,  i.e.,  ESF  plumes  that  have  stopped  rising,  is  exam¬ 
ined  using  the  NRL  SAMI3/ESF  three-dimensional  simula¬ 
tion  code.  We  find  that  fossil  bubbles,  plasma  density  de¬ 
pletions  associated  with  fossil  plumes,  can  persist  as  high- 
altitude  equatorial  depletions  even  while  being  “blown”  by 
zonal  winds.  Corresponding  airglow -proxy  images  of  fossil 
plumes,  plots  of  electron  density  versus  longitude  and  lati¬ 
tude  at  a  constant  altitude  of  288  km,  are  shown  to  partially 
“fill  in”  in  most  cases,  beginning  with  the  highest  altitude 
field  lines  within  the  plume.  Specifically,  field  lines  upon 
which  the  E  field  has  fallen  entirely  to  zero  are  affected  and 
only  the  low  altitude  (<600  km)  portion  if  each  field  line  fills 
in.  This  suggests  that  it  should  be  possible  to  observe  a  bub¬ 
ble  at  high  altitude  on  a  field  line  for  which  the  corresponding 
airglow  image  no  longer  shows  a  depletion.  In  all  cases  ESF 
plumes  stop  rising  when  the  flux-tube-integrated  ion  mass 
density  inside  the  upper  edge  of  the  bubble  is  equal  to  that  of 
the  nearby  background,  further  supporting  the  result  of  Krall 
et  al.  (2010b). 
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1  Introduction 

Equatorial  spread  F  (ESF),  the  result  of  a  Rayleigh-Taylor- 
like  instability  in  the  night  time  ionospheric  F-layer  (Haeren- 
del,  1974;  Ossakow,  1981;  Kelley,  1989;  Hysell,  2000),  has 
deleterious  affects  on  electromagnetic  signals.  For  example, 
ESF  was  discovered  because  it  causes  “diffuse  echoes  from 
the  /■’-region  of  the  ionosphere  received  continuously  at  night 
in  equatorial  regions  over  a  wide  range  of  wave -frequency” 


Correspondence  to:  J.  Krall 
(jonathan.krall@nrl.navy.mil) 


(Booker  and  Wells,  1938).  At  times,  the  strong  electron  den¬ 
sity  gradients  associated  with  ESF  degrade  global  commu¬ 
nication  and  navigation  systems  (de  La  Beaujardiere  et  al., 
2004;  Steenburgh  et  al.,  2008).  Recently,  the  use  of  sophisti¬ 
cated  three-dimensional  computer  simulation  models  (Huba 
et  al.,  2008;  Retterer,  2010)  has  increased  our  basic  under¬ 
standing  of  ESF.  Computer  simulations  have  provided  newly 
detailed  explanations  of  numerous  ESF-related  phenomena 
including  compositional  (Huba  et  al.,  2009b)  and  tempera¬ 
ture  signatures  (Huba  et  al.,  2009a)  of  ESF,  ESF  stabilization 
by  meridional  winds  (Krall  et  al.,  2009a;  Mamyama  et  al., 
2009;  Abdu  et  al.,  2006),  the  effect  of  zonal  winds  on  ESF 
plume  morphology  (Huba  et  al.,  2009c),  ESF  airglow  (Krall 
et  al.,  2009b),  and  ESF-related  “plasma  blobs”  (Krall  et  al., 
2010a). 

Radar  backscatter  observations  show  that  an  ESF  plume 
can  persist  for  hours  after  it  has  halted  its  upward  motion 
with  correspondingly  persistent  effects  on  co-located  elec¬ 
tromagnetic  signals  (Yokoyama  et  al.,  2007).  Plumes  that 
have  halted  their  upward  motion  are  called  “fossils”  (Argo 
and  Kelley,  1986;  Mendillo  et  al.,  1992;  Schunk  and  Sojka, 
1996;  Makela,  2006;  Sekar  et  al.,  2007). 

In  this  paper,  we  use  the  SAMI3/ESF  three-dimensional 
simulation  model  (Huba  et  al.,  2008)  to  study  the  properties 
and  signatures  of  fossil  ESF  plumes  and  the  effect  of  winds 
that  transport  these  phenomena  across  the  geomagnetic  field 
in  the  zonal  direction.  We  also  consider  the  case  in  which  a 
converging  zonal  wind  acts  to  halt  the  upward  rise  of  an  ESF 
bubble,  leading  to  a  fossil. 

This  study  is  a  follow-on  to  two  previous  studies.  In  one 
we  showed  that,  in  the  case  of  zero  neutral  winds,  an  ESF 
bubble  stops  rising  and  becomes  a  fossil  when  the  flux-tube- 
averaged  ion  mass  density  at  the  upper  edge  of  the  bubble 
rises  to  the  point  that  it  is  equal  to  the  corresponding  ambient 
value  (Krall  et  al.,  2010b).  We  further  showed  that  a  fossil 
bubble  continues  to  persist  and  evolve  for  hours  after  it  has 
stopped  rising.  In  another  (Krall  et  al.,  2009b),  we  found 
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that  a  converging  zonal  wind  can  lead  to  the  occurrence  of 
a  fossil  bubble  at  lower  heights  and  earlier  simulation  times 
than  would  otherwise  occur.  We  found  that  such  bubbles  are 
accompanied  by  a  unique  airglow  signature,  but  did  not  re¬ 
late  this  finding  to  other  modelling  or  observations  of  fossil 
plumes  or  bubbles.  The  aim  of  the  present  study  is  to  de¬ 
scribe  simulations  of  various  cases  in  which  fossil  plumes 
are  expected  to  occur  and  to  do  so  in  enough  detail  that  these 
findings  can  be  verified  through  simultaneous  in  situ,  radar, 
and/or  airglow  observations.  Here  we  are  particularly  inter¬ 
ested  in  the  effects  of  winds. 

We  find  that  ESF  plumes  that  have  stopped  rising  often 
contain  regions  with  persistent,  but  relatively  weak,  electric 
fields  and  that  these  regions  show  a  distinct  airglow  signa¬ 
ture.  In  particular,  field  lines  that  support  these  persistent  E 
fields  continue  to  show  airglow  depletions  even  while  other 
parts  of  the  same  ESF  airglow  depletion  are  “filling  in.”  We 
find  that  fossil  plumes  and  corresponding  plasma  density  de¬ 
pletions  are  transported  by  zonal  winds,  as  observed,  and  that 
the  density,  electric  field  and  airglow  signatures  found  in  the 
zero-wind  simulation  also  occur  for  wind-transported  ESF 
plumes.  We  show  that  a  plume  that  is  interrupted  by  converg¬ 
ing  zonal  winds,  resulting  in  a  fossil,  has  an  electric  field  that 
falls  everywhere  to  zero,  so  that  the  entirety  of  the  associ¬ 
ated  airglow  depletion  slowly  fills  in.  In  this  case  the  airglow 
depletion  might  appear  to  fade  out  entirely  while  a  plasma 
density  depletion  on  the  same  field  line  persists  at  higher  al¬ 
titudes.  Finally,  we  show  that  a  weakly  growing  plume  can 
attain  a  fossilized  state  where  the  E  field  is  everywhere  weak 
but  persistent,  resulting  in  a  weak  but  persistent  depletion  at 
the  airglow  layer. 

2  SAMI3  simulations 

The  NRL  3-D  ionosphere  code  SAMI3/ESF  (Huba  et  al., 
2008)  is  used  in  this  analysis.  SAMI3/ESF  is  a  version  of 
the  SAMI3  global  ionosphere  code  (Huba  et  al.,  2000,  2005) 
that  has  been  modified  for  the  study  of  ESF.  SAMI3  models 
the  plasma  and  chemical  evolution  of  seven  ion  species  (H+, 
He+,  N+,  0+,  N+,  NO+  and  07).  The  complete  ion  tem¬ 
perature  equation  is  solved  for  three  ion  species  (H+,  He+ 
and  0+)  as  is  the  electron  temperature  equation.  Ion  inertia 
is  included  in  the  ion  momentum  equation  for  motion  along 
the  geomagnetic  field  and  E  x  B  drifts  are  computed  to  ob¬ 
tain  motion  transverse  to  the  field.  Neutral  composition  and 
temperature  are  specified  using  the  empirical  NRLMSISE00 
model  (Picone  et  al.,  2002).  The  version  of  SAMI3  used 
here  and  in  our  other  recent  ESF  studies  (Huba  et  al.,  2008, 
2009a,b,c;  Krall  et  al.,  2009a,b,  2010a,b)  includes  a  per¬ 
pendicular  electric  field  E±  —  —VO  that  is  computed  self- 
consistently  to  determine  the  E  x  B  drifts. 

The  SAMI3/ESF  potential  equation  is  derived  from  cur¬ 
rent  conservation  (V  ■  /  =  0)  in  dipole  coordinates  (,v,  p,<f>) 
and  is  similar  to  that  derived  in  Haerendel  et  al.  (1992): 


3  30  3  30  dFpy 

dp  pp  dp  +  dtp  3 cp  3 (p  dp 


(1) 


where  F^g  =  / (rEshr'9/A)(B0/c)crHcgpds,  and  FpV  = 
f  rgsint? (Bo / c)o pVn^ds .  Here,  O  is  the  electrostatic  poten¬ 
tial,  gp  is  the  component  of  gravity  perpendicular  to  B,  op  is 
the  Pedersen  conductivity. 


a  He 


=E 
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is  the  ion  component  of  the  Hall  conductivity  divided  by 
the  ion  cyclotron  frequency,  Epp  —  f  (p A / bs)opds ,  Eprj,  = 
/ {\/pbs  A)opds,  A  =  (1  -(Geos2#)1/2,  G;  =eB/mjC,  vm  is 
the  ion-neutral  collision  frequency,  B  is  the  local  geomag¬ 
netic  field,  Bq  is  the  geomagnetic  field  at  the  equator,  9  is 
the  magnetic  colatitude,  bs  =  (r|/r3)  A,  rE  is  the  radius  of 
the  earth,  and  field-line  integrations  with  respect  to  coordi¬ 
nate  s  are  along  the  entire  field  line,  with  the  base  of  the 
field  lines  being  at  85  km.  Equation  (1)  has  been  simplified 
by  neglecting  Hall-conductance  terms  relative  to  Pedersen- 
conductance  terms  on  the  LHS  (e.g.,  'Em  —  oh  —  0),  neglect¬ 
ing  the  dFpg/dp  term  relative  to  the  dF^g/dcp  term  on  the 
RHS,  and  by  setting  the  meridional  wind  to  zero.  Examples 
of  the  potential  equation  without  these  simplifications  appear 
in  Krall  et  al.  (2009b)  and,  in  the  context  of  the  SAMI3  code, 
Huba  et  al.  (2010). 

The  3-D  simulation  model  is  initialized  using  results 
from  the  two-dimensional  SAME  code.  SAME  is  run  for 
48  h  using  the  following  geophysical  conditions:  F10.7  =  75, 
F10.7A  =  75,  Ap  =  4,  and  day-of-year  80  (e.g.,  21  March 
2002).  The  simulation  is  centered  at  geographic  longitude 
0°  so  universal  time  and  local  time  are  the  same.  The  E  x  B 
drift  in  SAME  is  prescribed  by  the  Fejer/Scherliess  model 
(Scherliess  and  Fejer,  1999).  The  plasma  parameters  (den¬ 
sity,  temperature,  and  velocity)  at  time  19:20  UT  (of  the  sec¬ 
ond  day)  are  used  to  initialize  the  3-D  model  at  each  mag¬ 
netic  longitudinal  plane. 

The  3-D  model  uses  a  non-uniform,  non-orthogonal  nu¬ 
merical  grid  that  follows  the  geomagnetic  field.  The  simu¬ 
lation  region  includes  field  lines  with  magnetic  apex  heights 
from  90  km  to  2400km,  and  longitudes  from  —2°  to  2°  (lon¬ 
gitudinal  width  of  4°  or  ~  460  km).  The  plasma  is  modelled 
from  hemisphere  to  hemisphere  up  to  ±31°  magnetic  lati¬ 
tude.  Except  where  noted,  the  grid  is  (nz,«f,«i)  =  (101,  202, 
96),  where  nz  is  the  number  of  grid  points  along  the  magnetic 
field,  «f  is  the  number  of  “field  lines,”  and  n\  is  the  number 
of  longitudes.  This  grid  has  a  resolution  of  ~6  km  x  5  km  in 
altitude  and  longitude  in  the  magnetic  equatorial  plane.  The 
numerical  grid  is  periodic  in  longitude.  In  essence  we  are 
simulating  a  narrow  “wedge”  of  the  ionosphere  in  the  post¬ 
sunset  period.  Unless  noted,  the  ESF  instability  in  each  simu¬ 
lation  was  initiated  by  imposing  a  Gaussian-like  perturbation 
in  the  plasma  density  at  the  initial  time  (19:20  UT):  a  peak 
density  perturbation  of  15%  centered  at  0°  longitude  with  a 
half-width  of  0.25°,  and  at  an  altitude  z  =  300  km  with  a  half¬ 
width  of  50  km.  This  seed  perturbation  extends  along  the 


Ann.  Geophys.,  28,  2059-2069,  2010 


www.  ann-geophys  .net/28/2059/20 1 0/ 


J.  Krall  et  al.:  ESF  fossil  plumes 


2061 


E 


4) 

"D 

3 


UT  01:17 


UT  02:34 


-20  -10  0  10  20  -1  0  1  -2-1  0  1  2  -10  1 
Latitude  Longitude  Longitude  Longitude 


Fig.  1.  Left:  Contours  of  logins  versus  height  and  latitude  at  longitude  0  at  various  times.  Right  three  columns:  Contours  of  the  log  of 
the  flux-tube  averaged  ion  mass  density  (logjo(p/™o+))-  t^le  l°g  °f  the  Pedersen  conductance  (logio(Spp))  and  the  electrostatic  potential, 


plotted  in  the  equatorial  plane.  Number  densities  for  ne,  and  (p/mQ+),  are  in  units  of  cm  3  and  correspond  to  the  color-bar  scale.  In  the 
left  hand  panels,  0+  ion  velocity  vectors  indicate  the  direction  of  flow  away  from  the  dot.  Here  there  are  no  winds. 


entire  flux  tube.  For  this  study,  the  vertical  and  meridional 
neutral  winds  are  set  to  zero  and  the  zonal  neutral  wind  is  set 
either  to  zero,  to  a  constant  100  m  s_  1  or  to  a  time-dependent 
converging  wind  pattern  with  peak  velocities  ±  10  m  s“ 1 . 


3  The  development  of  a  fossil  bubble 

Figure  1  shows  an  example  of  a  simple  ESF  bubble  viewed 
in  a  constant  longitude  plane  (left  panels).  Here  the  zonal 
wind  is  set  to  zero.  The  smaller  right-hand  panels  in  this  fig¬ 
ure  show  contours  of  the  log  of  the  flux-tube  averaged  mass 
density  normalized  to  the  0+  mass  \ogw(p/mo+ ),  the  log 
of  the  Pedersen  conductance  log (T,pp),  and  the  electrostatic 
potential  <J>  plotted  in  the  equatorial  plane.  The  average  mass 
density,  which  is  reduced  within  an  active,  upwardly  drifting 
bubble,  indicates  the  shape  and  extent  of  the  bubble  that  is 


shown  as  contours  of  log10(«e)  in  the  left  hand  panels.  The 
left-hand  panel  shows  the  motion  of  0+  ions  (both  E  x  B 
drifts  and  field-aligned  flows),  where  velocity  vectors  indi¬ 
cate  the  direction  of  flow  away  from  the  dot  in  each  case.  In 
this  simulation  there  are  no  winds  ( Vn(p  —  0). 

This  simulation  was  described  in  Krall  et  al.  (2010a), 
where  it  was  shown  that  the  upward  motion  of  the  bubble 
stops  when  (p)  at  the  upward  edge  of  the  bubble  is  the  same 
as  that  of  the  nearby  background  ionosphere  at  the  same 
height.  Here  we  describe  this  result  in  further  detail.  We  first 
note  the  sequence  of  events  shown  in  Fig.  1.  In  the  upper 
row  of  contour  plots  (00:00  UT),  the  bubble  is  active.  This  is 
indicated  by  ion  flow  vectors  that  extend  outward  across  field 
lines  (the  dots  in  each  left-hand  panel  are  arrayed  along  geo¬ 
magnetic  field  lines),  by  (p)  values  within  the  bubble  that  are 
lower  than  background  values,  and  by  a  strong  zonal  gradient 
in  the  electrostatic  potential.  As  noted  in  Huba  et  al.  (2009b) 
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Fig.  2.  Electron  density  versus  longitude  and  latitude  at  a  constant  altitude  of  288  km  at  various  times.  The  first  three  time  values  correspond 
to  those  shown  in  Fig.  1.  The  latitude  range  lies  in  northern  half  of  the  simulation. 


and  Krall  et  al.  (2010a),  a  “super  fountain  effect”  acts  within 
an  ESF  bubble,  producing  density  crests  at  altitude  600  km 
and  latitude  ±16°  that  are  analogous  to  the  background  “Ap¬ 
pleton  anomaly”  ionization  crests.  In  the  middle  row  of  con¬ 
tour  plots  (01:17  UT),  the  bubble  has  nearly  halted  its  upward 
motion.  This  is  indicated  by  ion  flow  vectors  near  the  top  of 
the  bubble  that  are  almost  entirely  parallel  to  the  geomag¬ 
netic  field  (on  lower  altitude  field  lines,  upward  E  x  B  drifts 
are  more  prevalent),  by  (p)  at  the  top  of  the  bubble  that  is 
close  to  background  values,  and  by  a  nearly-zero  zonal  gra¬ 
dient  in  the  electrostatic  potential  at  the  height  of  the  bub¬ 
ble  apex.  Here  the  density  crests  within  the  bubble  (altitude 
500km,  latitude  ±18°)  have  fallen  to  lower  altitude  and  have 
grown  stronger  as  plasma  drains  from  the  bubble  apex. 

In  the  lower  row  of  Fig.  1  (02:34  UT),  the  bubble  has 
reached  a  largely-quiescent  fossil  state,  with  minimal  plasma 
motion  in  the  upper  part  of  the  bubble,  where  a  density  de¬ 
pletion  persists.  On  these  outer  field  lines,  which  we  refer 
to  as  the  “head”  of  the  bubble,  the  plasma  density  at  alti¬ 
tudes  of  300-700  km  is  now  similar  to  that  of  the  background 
ionosphere.  As  a  result,  the  Pedersen  conductance  Y.pp  is 
also  nearly  equal  to  that  of  the  background  ionosphere.  On 
lower  altitude  field  lines  within  the  bubble,  the  “stalk”  of  the 
bubble,  we  find  that  a  weak  electric  field  persists  along  with 
low  values  of  (p)  and  T,pp,  relative  to  the  background.  We 


note  that  the  7  h  of  simulation  time  that  have  passed  in  order 
for  the  bubble  to  obtain  a  fossil  state  is  probably  an  over¬ 
estimate.  Simulations  of  bubbles  in  a  global  context  by  Huba 
and  Joyce  (2010)  show  that  perturbations  that  are  imposed  as 
early  as  17:00  local  time  can  lead  to  growing  ESF  bubbles. 

Further  results  from  this  same  simulation  are  shown  in 
Fig.  2.  As  a  proxy  for  the  airglow  images  that  have  been 
used  in  the  past  to  study  ESF  depletions  (see,  e.g.,  Makela, 
2006),  we  show  contours  of  electron  density  versus  longi¬ 
tude  and  latitude  at  a  constant  altitude  of  288  km  at  various 
times,  the  first  three  of  which  correspond  to  those  of  Fig.  1. 
The  altitude  288  km  is  used  because  it  lies  in  the  altitude 
range,  250-300  km,  that  is  typically  the  source  of  630.0  nm 
airglow.  The  latitude  range  of  Fig.  2  is  in  the  northern  half 
of  the  simulation.  As  the  simulated  ESF  depletion  increases 
in  height,  it  extends  northward  (and  southward)  in  latitude, 
moving  further  into  the  field  of  view  of  Fig.  2.  As  the  bubble 
obtains  a  fossil  state  it  first  stops  moving  outward  in  latitude 
(01:17UT).  The  highest  latitude  areas  of  depletion  then  be¬ 
gin  to  fill  in  (02:34  UT)  as  the  plasma  density  within  the  bub¬ 
ble  falls  further  downward  (see  Fig.  1,  lower  left  panel).  At 
the  end  of  the  simulation  (03: 18  UT)  the  upper  part  (head)  of 
the  airglow  depletion  has  filled  in  while  the  airglow  signature 
of  the  stalk  remains  dark. 
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Fig.  3.  Left:  Contours  of  logj gne  versus  height  and  latitude  at  various  times.  Right  three  columns:  Contours  of  logiQUe,  the  log  of  the 
flux-tube  averaged  ion  mass  density  (log1o(/o/mQ+))  and  the  log  of  the  Pedersen  conductance  (logio(Ep/>)),  plotted  in  the  equatorial  plane. 
Number  densities  for  ne,  and  (p/m0+),  are  in  units  of  cm-3  and  correspond  to  the  colour-bar  scale.  In  the  left  hand  panels,  0+  ion  velocity 
vectors  indicate  the  direction  of  flow  away  from  the  dot.  This  case  includes  a  constant,  eastward  zonal  wind  of  100  nr  s-  1 . 


Our  basic  findings  are  that  the  airglow  signature  can  decay 
even  while  the  bubble  persists  and  that  a  persistent  airglow 
signature  is  associated  with  persistent,  if  weak,  E  fields.  The 
collapse  of  the  electric  field  in  the  head  of  the  bubble  at  late 
times  is  associated  with  a  “filling  in”  of  both  the  Pedersen 
conductance  and  the  simulated  airglow  depletion. 


4  The  effect  of  a  zonal  wind  on  a  fossil  bubble 

Figures  3  and  4  show  a  repeat  of  the  simulation  of  Figs.  1 
and  2  but  with  a  constant  eastward  100ms-1  zonal  neutral 
wind  imposed.  As  described  by  Zalesak  et  al.  (1982),  zonal 
winds  act  on  the  density  gradients  of  an  ESF  bubble  so  as 
to  transport  the  bubble  in  the  direction  of  the  wind,  albeit 
with  a  velocity  below  that  of  the  wind  velocity.  In  order  to 
accurately  capture  this  transport,  the  number  of  zonal  grid 
points  was  increased  from  96  to  192  for  this  simulation. 


Shown  in  Fig.  3  is  the  bubble  viewed  in  a  constant  longi¬ 
tude  plane  (left  panels)  as  contours  of  log10(«e)  at  various 
times.  As  in  Fig.  1,  0+  motion  is  shown  by  velocity  vectors 
that  indicate  the  direction  of  flow  away  from  the  dot  in  each 
case.  The  smaller  right-hand  panels  in  this  figure  show  con¬ 
tours  of  the  log  of  the  electron  density  in  the  equatorial  plane 
log|0nf,,  the  flux-tube  averaged  mass  density  normalized  to 
the  0+  mass  \og[0(p/mQ+),  and  the  log  of  the  Pedersen  con¬ 
ductance  logfEpp).  The  location  in  longitude  for  each  left- 
hand  panel  is  indicated  by  a  vertical  line  in  the  corresponding 
plot  of  log10He  versus  longitude  and  altitude.  Figure  3  shows 
bubbles  that  are  tilted  towards  the  western  direction  in  com¬ 
parison  to  the  zero- wind  result  of  Fig.  1,  which  shows  no  tilt. 
Without  a  realistic  wind  pattern,  as  was  used  in  Huba  et  al. 
(2009c),  commonly-observed  “C”  shaped  plumes  are  not  ob¬ 
tained. 

Further  results  from  this  same  simulation  are  shown  in 
Fig.  4,  where  contours  of  electron  density  versus  longitude 
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Fig.  4.  Electron  density  versus  longitude  and  latitude  at  a  constant  altitude  of  288  km  at  various  times.  The  first  three  time  values  correspond 
to  those  shown  in  Fig.  3. 


and  latitude  are  plotted  at  a  constant  altitude  of  288  km  at 
various  times,  the  first  three  of  which  correspond  to  those  of 
Fig.  3.  Time  values  for  the  three  rows  of  Fig.  3  and  the  first 
three  panels  of  Fig.  4  were  chosen  to  show  configurations 
similar  to  those  of  Fig.  1 :  an  active  bubble,  a  bubble  that  has 
nearly  halted,  and  a  fossil  bubble  that  has  evolved  for  some 
time. 

The  key  result  of  this  simulation  is  that  the  findings  in  the 
simple  zero-wind  simulation  also  hold  for  the  more  realistic 
case  where  a  zonal  wind  transports  the  fossil  bubble  even 
as  vertical  and  field-aligned  ion  motion  within  the  bubble 
are  decaying,  first  in  the  head  of  the  plume  and  then  in  the 
stalk.  As  expected,  the  presence  of  the  zonal  wind  results 
in  a  bubble  that  does  not  rise  as  quickly  as  in  the  zero-wind 
case  (Zalesak  et  al.,  1982).  As  in  the  zero-wind  simulation, 
the  airglow  signature  of  the  head  fades  as,  simultaneously, 
the  Pedersen  conductance  increases  and  the  plasma  drifts  fall 
virtually  to  zero. 

As  is  seen  in  radar  observations,  the  fossil  density  deple¬ 
tion  above  the  F-layer  persists  even  as  it  is  distorted  by  its 
interaction  with  the  zonal  wind.  Because  the  simulation  does 
not  use  a  high-order,  non-diffusive  transport  scheme,  neither 
the  active  (e.g.,  21:30  UT)  nor  the  fossil  depletion  show  the 


bifurcations  that  are  evident  in  such  simulations  (see,  e.g., 
Huba  and  Joyce,  2007). 

5  The  effect  of  a  converging  zonal  wind 

As  shown  in  Figs.  5  and  6,  we  repeat  once  again  the  simula¬ 
tion  of  Figs.  1  and  2  but  with  a  very  mild  (10  m  s_1)  converg¬ 
ing  zonal  wind  imposed,  beginning  at  22:00  UT  and  reach¬ 
ing  full  force  at  23:00  UT.  The  wind  is  zero  at  the  edges  of 
the  simulation  (±2°  longitude)  and  ramps  up  to  full  speed 
(±10m  s'1)  at  the  center  of  the  simulation.  In  other  words, 
the  wind  is  eastward  for  longitudes  <  0  and  westward  for  lon¬ 
gitudes  >  0  (the  wind  pattern  does  not  vary  with  height  and 
corresponding  meridional  winds  are  neglected).  This  wind 
pattern  is  simular  to  that  used  in  Krall  et  al.  (2009b)  and,  as 
in  that  study,  the  wind  has  the  effect  of  halting  the  upward 
motion  of  the  ESF  bubble  at  an  earlier  time  and  a  lower  alti¬ 
tude  than  in  an  otherwise-identical  simulation  with  no  winds. 

Shown  in  Fig.  5  is  the  bubble  viewed  in  a  constant  longi¬ 
tude  plane  (left  panels)  as  contours  of  \ogl0(ne).  The  smaller 
right-hand  panels  show  contours  of  the  flux-tube  averaged 
mass  density  normalized  to  the  0+  mass  log10(p/;n0+), 
the  log  of  the  Pedersen  conductance  log(Ew),  and  the 
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Fig.  5.  Left:  Contours  of  log iQtie  versus  height  and  latitude  at  longitude  0  at  various  times.  Right  three  columns:  Contours  of  the  log  of 
the  flux-tube  averaged  ion  mass  density  (logj0(p/mQ+)),  the  log  of  the  Pedersen  conductance  (logiotEpp))  and  the  electrostatic  potential, 
plotted  in  the  equatorial  plane.  Number  densities  for  ne,  and  ( p/irq+ ),  are  in  units  of  cm_J  and  correspond  to  the  colour-bar  scale.  In  the 
left  hand  panels,  0+  ion  velocity  vectors  indicate  the  direction  of  flow  away  from  the  dot.  This  case  includes  a  converging  zonal  wind  of 
±10  ms-1. 


electrostatic  potential  <J>  plotted  in  the  equatorial  plane.  As 
in  Figs.  1  and  3,  0+  motion  is  plotted  via  velocity  vectors 
that  indicate  the  direction  of  flow  away  from  the  dot  in  each 
case. 

Airglow-proxy  results  from  this  simulation  are  shown  in 
Fig.  6,  where  contours  of  electron  density  versus  longitude 
and  latitude  are  plotted  at  a  constant  altitude  of  288  km  at 
various  times,  the  first  three  of  which  correspond  to  those  of 
Fig.  5.  Once  again,  time  values  for  the  three  rows  of  Fig.  5 
and  the  first  three  panels  of  Fig.  6  were  chosen  to  show  spe¬ 
cific  states  of  the  evolving  bubble:  an  active  bubble,  a  bubble 
that  has  nearly  halted,  and  a  fossil  bubble  that  has  evolved 
for  some  time. 

The  key  difference  between  these  results  and  those  of 
Figs.  1  and  2  is  that  the  converging  wind  has  had  the  effect 
of  reducing  and  then  stopping  upward  E  x  B  drifts  through¬ 


out  the  bubble;  in  the  other  cases,  this  effect  is  limited  to  the 
head  of  the  bubble.  Figure  6  notably  differs  from  Figs.  2  and 
4  in  that  the  airglow  signature  of  the  entire  depletion  fades 
instead  of  only  the  upper  part.  As  in  the  other  cases  above, 
the  plasma  density  depletion  at  high  altitudes  persists  as  the 
bubble  becomes  a  fossil. 

6  Fossils  resulting  from  weakly-growing  bubbles 

Finally,  we  repeat  the  simulation  of  Figs.  1  and  2  but  with  the 
instability  seeded  via  random  fluctuations  in  the  0+  density, 
temperature,  and  velocity.  As  in  Krall  et  al.  (2010a)  these 
initial  fluctuations  are  based  on  the  observed  amplitudes  of 
traveling  ionospheric  disturbances  as  reported  in  Hocke  and 
Schlegel  (1996,  see  Fig.  17  therein).  This  simulation  was 
previously  described  in  Krall  et  al.  (2010a)  as  an  example  of 
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Fig.  6.  Electron  density  versus  longitude  and  latitude  at  a  constant  altitude  of  288  km  at  various  times.  The  first  three  time  values  correspond 
to  those  shown  in  Fig.  5. 


the  unusual  situation  where  the  ESF  “super  fountain”  (Huba 
et  al.,  2009b)  is  not  very  strong,  leading  to  “bubbles”  that 
have  mild  plasma  density  enhancements  (“plasma  blobs,” 
see,  e.g.,  Krall  et  al.,  2010a,  and  references  therein)  instead 
of  strong  depletions  at  high-altitude. 

Results  are  shown  in  Figs.  7  and  8  where  a  selected 
bubble  is  viewed  in  a  constant  longitude  plane  as  con¬ 
tours  of  \oglQ(ne)  (Fig.  7,  left-hand  panels),  as  contours  of 
the  flux-tube  averaged  mass  density  normalized  to  the  0+ 
mass  log10(p/fl?0+}.  the  log  of  the  Pedersen  conductance 
log (Epp),  and  the  electrostatic  potential  <t>  plotted  in  the 
equatorial  plane  (Fig.  7,  small  panels),  and  as  contours  of 
electron  density  versus  longitude  and  latitude  are  plotted  at  a 
constant  altitude  of  288  km  (Fig.  8).  As  in  Fig.  1,  0+  motion 
is  shown  in  Fig.  7  as  velocity  vectors  that  indicate  the  direc¬ 
tion  of  flow  away  from  the  dot  in  each  case.  As  in  Fig.  3, 
the  location  in  longitude  for  each  left-hand  panel  of  Fig.  7,  is 
indicated  by  a  vertical  line  in  one  of  the  corresponding  plots 
versus  longitude  and  altitude.  As  in  the  previous  cases,  time 
values  for  the  three  rows  of  Fig.  7  and  the  first  three  panels 
of  Fig.  8  were  chosen  to  show  specific  states  of  the  evolving 
bubble:  an  active  bubble,  a  bubble  that  has  nearly  halted,  and 
a  fossil  bubble  that  has  evolved  for  some  time. 


As  in  all  other  cases,  the  upward  motion  of  the  ESF  plume 
halts  when  (p)  values  within  the  apex  of  the  plume  are  equal 
to  those  in  the  nearby  background.  Consistent  with  Chen 
et  al.  (1983),  an  ESF  plume  in  a  system  with  multiple  bub¬ 
bles  rises  more  slowly  than  in  the  single-bubble  case  (Chen 
et  al.,  1983,  also  showed  that  the  fields  in  a  multi-bubble  sys¬ 
tem  tend  to  distort  the  bubbles,  as  is  seen  here).  As  a  result  of 
this  weak  ESF  growth,  the  super  fountain  effect  is  not  strong 
(here  the  scale  of  the  plotted  vectors  is  only  100  ms-1)  and 
there  is  no  plasma  density  depletion  at  the  apex  of  the  fos¬ 
sil  plume.  Because  of  this  lack  of  a  high-altitude  depletion, 
and  because  we  see  this  behaviour  only  in  simulations  with 
both  very  low  F10J  and  F10JA  indices  and  with  multiple 
plumes,  we  believe  this  case  to  be  somewhat  unusual.  Nev¬ 
ertheless,  this  case  is  interesting  in  that  both  the  E  field  and 
the  full  height  of  the  airglow  depletion  persist  through  to  the 
near-dawn  end  of  the  simulation.  This  supports  the  hypothe¬ 
sis,  suggested  by  the  previous  three  simulations,  that  airglow 
signatures  are  associated  with  field  lines  where  the  E  field  is 
nonzero. 
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Fig.  7.  Left:  Contours  of  log]gne  versus  height  and  latitude  at  longitude  1.6°  at  various  times.  Right  three  columns:  Contours  of  the  log  of 
the  flux-tube  averaged  ion  mass  density  (logjo(p/mQ+)),  the  log  of  the  Pedersen  conductance  (logio(Spp))  and  the  electrostatic  potential. 


plotted  in  the  equatorial  plane.  Number  densities  for  ne,  and  ( p/niQ+ ),  are  in  units  of  cm--’  and  correspond  to  the  colour-bar  scale.  In  the 
left  hand  panels,  0+  ion  velocity  vectors  indicate  the  direction  of  flow  away  from  the  dot.  Here  there  are  no  winds  and  the  instability  is 
seeded  by  imposing  random  fluctuations. 


7  Conclusions 

We  have  used  the  NRL  SAMI3/ESF  simulation  code  to  ex¬ 
amine  the  behaviour  of  fossil  bubbles  in  simulations  in  which 
we  varied  either  the  zonal  neutral  winds  or  the  initial  con¬ 
dition  for  the  instability.  These  support  the  previous  re¬ 
sult  of  Krall  et  al.  (2010b):  in  all  simulations,  ESF  plumes 
stopped  rising  when  the  flux-tube-integrated  ion  mass  den¬ 
sity  inside  the  upper  edge  of  the  bubble  was  equal  to  that  of 
the  nearby  background.  We  further  showed  that  fossils  can 
persist  as  high-altitude  plasma  density  depletions  even  while 
being  “blown”  by  a  zonal  wind.  Corresponding  airglow  im¬ 
ages  of  these  bubbles,  typically  observed  as  light  emitted  at 
altitudes  of  250-300  km  (corresponding  to  a  wavelength  of 
630.0  nm),  are  shown  to  fill  in  as  the  E  field  within  the  bub¬ 
ble  decays  to  zero  and  plasma  density  “falls”  down  field  lines 


to  low  altitudes.  In  particular,  field  lines  that  fill  in  are  those 
upon  which  the  E  field  has  fallen  entirely  to  zero  while  field 
lines  that  do  not  fill  in  support  a  persistent,  if  weak,  E  field. 
Typically,  only  the  higher  latitude  portion  of  the  airglow  de¬ 
pletion  signature  fills  in.  However,  a  plume  that  has  had  its 
upward  motion  interrupted  by  a  converging  zonal  wind  pat¬ 
tern  can  have  an  airglow  signature  that  fills  in  entirely.  In 
many  cases,  it  is  possible  to  have  a  fossil  plasma  density  de¬ 
pletion  at  high  altitude  (e.g.,  above  the  F-layer)  on  a  field  line 
for  which  the  corresponding  airglow  image  shows  no  deple¬ 
tion. 

In  these  simulations,  completely  decayed  E  fields  are  as¬ 
sociated  with  field  lines  where  the  plasma  has  fallen  back 
to  low  altitudes,  filling  in  the  airglow  images  and  increas¬ 
ing  the  Pedersen  conductance.  In  related  work.  Park  et  al. 
(2010)  has  considered  the  reduction  in  the  E  field  by  a 
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Fig.  8.  Electron  density  versus  longitude  and  latitude  at  a  constant  altitude  of  288  km  at  various  times.  The  first  three  time  values  correspond 
to  those  shown  in  Fig.  7. 


high-conductance  region  in  order  to  explain  observations  of 
“plasma  blobs.”  In  that  work,  there  is  evidence  that  field- 
aligned  currents  are  partially  closed  across  blobs,  which  are 
high  density  regions  within  an  ESF  plume  (see  Krall  et  al., 
2010a,  and  references  within).  This  leads  us  to  consider  the 
idea  that  a  localized  region  within  an  ESF  plume  with  a  high 
value  of  the  Pedersen  conductivity  can  “short  out”  the  E  field 
on  coincident  field  lines  within  the  plume,  affecting  the  dy¬ 
namics  of  the  plume.  However,  based  on  the  sequence  of 
events  in  each  simulation,  where  the  bubble  stops  rising  be¬ 
cause  of  conditions  at  the  outward  leading  edge  of  the  bubble 
even  while  E  fields  within  the  bubble  are  quite  large,  we  con¬ 
clude  that  any  “shorting  out”  of  the  field  by  an  increased  Ped¬ 
ersen  conductance  is  a  secondary  effect  that  occurs  after  the 
bubble  has  stopped  rising  and  has  reached  a  slowly-evolving 
fossil  state. 
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